A non-monotonic constitutive model is not necessary to obtain shear banding 
phenomena in entangled polymer solutions 
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In 1975 Doi and Edwards predicted that entangled polymer melts and solutions can have a con- 
stitutive instability, signified by a decreasing stress for shear rates greater than the inverse of the 
reptation time. Experiments did not support this, and more sophisticated theories incorporated 
Marrucci's idea (1996) of removing constraints by advection; this produced a monotonically increas- 
ing stress and thus stable constitutive behavior. Recent experiments have suggested that entangled 
polymer solutions may possess a constitutive instability after all, and have led some workers to ques- 
tion the validity of existing constitutive models. In this Letter we use a simple modern constitutive 
model for entangled polymers, the Rolie-Poly model with an added solvent viscosity, and show that 
(1) instability and shear banding is captured within this simple class of models; (2) shear banding 
phenomena is observable for weakly stable fluids in flow geometries that impose a sufficiently inho- 
mogeneous total shear stress; (3) transient phenomena can possess inhomogeneities that resemble 
shear banding, even for weakly stable fluids. Many of these results are model-independent. 



PACS numbers: 83.80.Rs, 83.10.Kn, 83.60.Wc, 83.10.Gr 

Much of the rheology of entangled polymers solutions 
and melts is captured by the molecular theory of Doi 
and Edwards (DE) 0, who argued that polymers relax 
by curvilinear diffusion (reptation) within a tube of the 
surrounding polymers. The DE model has a local maxi- 
mum in the constitutive relation (the total shear stress as 
a function of shear rate for homogeneous flows). The re- 
sulting non-monotonic relation {e.g. the dashed curves in 
Fig. \2§ leads to an instability that for many years was not 
observed in experiments , but nonetheless attracted at- 
tention [HQ. This stress maximum is predicted to be less 
pronounced or absent if the convected constraint release 
(CCR) of entanglements due to flow 0, @, 0] is incorpo- 
rated. A sudden release of a constraint can relax both 
the orientation and conformation of a stretched polymer, 
which increases the stress and, for sufficiently frequent 
events, eliminates the instability. The CCR mechanism 
also leads to neutron scattering predictions that agree 
with experiment 0]. Similar physics applies to solutions 
of breakable wormlike micelles, in which the instability is 
well documented experimentally and leads to shear band- 
ing, in which a fast-flowing oriented state coexists with a 
more disordered and viscous state along a stress plateau 
There, CCR is less pronounced because of breakage 
and fails to ameliorate a constitutive instability 

Recently, Wang, Hu and co-workers studied entangled 
solutions of a high molecular weight (HMW) polymer 



11 hj, or DNA solutions 



in its own oligomer [ltji 11 . 
14| . finding a number of results that may be consistent 



with instability and shear banding after all. In controlled 
shear rate mode a weakly increasing stress plateau of 
three decades in shear rate was found, whereas in con- 
trolled shear stress mode the sheared solution experi- 



enced a jump in the shear rate, together with spatially 
inhomogeneous birefringence [lflj. Local velocimetry re- 
vealed spatially inhomogeneous velocity profiles in both 
the transient and the steady state [ijj regimes, while 
large amplitude oscillatory shear flow (LAOS) experi- 
ments showed an inhomogeneous banding-like shear rate 
profile at finite frequencies [l5| . Similar behavior was 
observed in a sliding plate shear cell in monodisperse so- 
lutions [lil ]. Relaxation after a step strain induced a 
highly inho mog eneous velocity field with negative local 



shear rates |13l |. Hu et al. found similar inhomogeneous 



flow behavior and possible signatures of shear banding in 
polymer solutions, and wormlike micelle solutions at con- 
centrations where severe shear thinning, but not banding, 
might be expected flrjj |. 

Wang et al. could not reconcile their results with exist- 
ing theory, and proposed that the instability is a yield like 
effect due to an unbalanced "entropic retraction force" 
17| . Here we show that much of the phenomenology of 
these experiments is consistent with the predictions of 
tube models with CCR, perhaps as anticipated in the 
original theory without introducing new physics. 

Model — We approximate the total stress T as separat- 
ing into fast Newtonian (or solvent) degrees of freedom 
and a slow viscoelastic component X (HMW polymer): 



T = -pi + 2nB + GS, 



(1) 



where I is the identity tensor, D = \ [Vv + (Vv) T ] , p 
is the isotropic pressure determined by incompressibility 
(V • v = 0), and rj is the solvent viscosity, for which we 
use the dimensionless quantity e = ij/(GTd). Here, Td 
is the reptation time. We are interested in the creeping 
flow (low Reynolds number) regime, in which V • T = 0. 
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FIG. 1: a) Parameters (/3,e) where the constitutive curve of 
the Rolie-Poly model is non-monotonic (shaded), b) Section 
through Couette rheometer showing the flow field and total 
stress components, T. 



For the dynamics of X we use the Rolie-Poly (RP) 
model [||, a simplified tube model that incorporates CCR 
0, where S(r, i) obeys 



(d t + v ■ V)S + (Vv) • S + £ • (Vv 
2 



\ T + — S 

Td 



2D — —(1—A)[I + S(l + f3A)] + DV S, (2) 

,4 = (l+trS/3) -1 / 2 and the Rouse time tr governs chain 
stretch. The stress "diffusion" term X>V 2 X describes the 
response to an inhomogencous viscoclastic stress; while 
not in the original RP model, it can arise due to diffusion 
or finite persistence length [la, il. We sp ecify Neu- 
mann boundary conditions (VS = 0) 0, Hlj. From ex- 
perimental values of the plateau modulus G ~ 6 x 10 2 Pa, 
reptation time r d ~ 20 s, and solvent viscosity r\ ~ 1 Pa s 
[10j . we use e = 10~ 5 . Here we use Td/r^ ~ 10 3 , which is 
consistent with the length of the stress-shear rate plateau 
reported in [T3| ■ The parameter (3 controls the efficiency 
of CCR, and is difficult find a precise value of experimen- 
tally. Ref. Q chose (3 — 1 to fit steady state data in poly- 
mer melts, and used multiple modes with f3 = 0.5 to fit 
experimental transient data. Here we tune between two 
qualitatively diffcrrent types of constitutive curve; either 
a non-monotonic (0.65) or a monotonic (0.728) constitu- 
tive curve with a broad plateau (Figs. [TT2|) . 

Eq. j2|) was solved in one spatial dimension using the 
Crank-Nicolson algorithm [la], for unidirectional Cou- 
ette flow v(r,t)6 between cylinders of radii i?i and i?2 
parameterized by q = In R2/ R\ . In this geometry the 
total shear stress T r g ~ 1/r 2 , so that the stress differ- 
ence across the flow cell is A\n.T r g = 2q. Cone and 
plate flow with cone angles of 9 = (4°, 1°) has been re- 
ported 11 , IH , so we use consistent values of stress differ- 
ence corresponding to q ~ AR/R =(2x 10 -3 , 2 x 10~ 4 ) 
[2(ij ]. Stresses are measured in units of G, shear rates in 
units of rj" 1 , and velocities in units of qRi/r « AR/tj 
for small q. To plot numerical data we use T, the di- 
mcnsionless specific torque (per height per radian) on 



the inner cylinder. The diffusion constant used was 
VT d /(R iq ) 2 = 4 x 10- 4 . 

Flow Curves — To calculate the steady state flow curves 
a step shear rate was applied from rest and evolved for 
500rd with time step 10~ 5 T d , after which subsequent 
shear rate steps and time evolutions were applied to scan 
up and down in shear rate (Fig. [5]). For non-monotonic 
constitutive curves (f3 = 0.65) shear banding always oc- 
curs, with hysteresis and a stress "plateau" . For the 
monotonic case (/3 = 0.728) shear banding could be in- 
ferred in the more highly curved geometry with the larger 
stress difference (larger q) , since the flow curve no longer 
follows the constitutive curve; but without hysteresis. 
Crudely, a monotonic flow curve exhibits banding-like 
flows when most of the shear rate in the gap occurs over 
a small range of stresses, i.e. the slope of the plateau 
must be much smaller than the apparent slope specified 
by the flow geometry: 



dr 



rORO-ro^) 



c.c. 



A7 



1, 



(3) 



where "C.C." denotes the flat portion of the constitutive 
curve (dashed in Fig. [2]) and "g" refers to the range of 
torques and shear rates specified by the flow geometry. 
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FIG. 2: Flow curves (solid black) and constitutive curves 
(dashed red) for a range of stress gradients q and CCR values 
j3. The shaded area shows the size of the torque difference 
Alog lo r = log 10 e 2<J . Circles indicate applied shear rates at 
which the transient state responses are shown in Fig. [3] 



The steady state velocity profiles are shown in Fig. [3] 
as solid (red) lines. The non-monotonic flow curves 
((3 = 0.65) lead to a pronounced kink in the velocity 
profile, a signature of shear banding. The monotonic 
case does not shear band in the flatter geometry (small 
q) , but for a more curved geometry (larger q) more shear 
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rates are accessible and the resulting smooth velocity pro- 
file could easily be interpreted as banding [ill ]; certainly 
the constitutive curve is not followed (Fig. [J). Similar 
smooth profiles were reported in [111 [l6j| , in a flow geom- 
etry with q ~ 0.004 — 0.02. A slightly increasing stress 
plateau over several decades in shear rates (as in [Io| |) 
would thus lead to apparently banding (inhomogenous) 
flow in geometries with very low stress gradients; but a 
linear steady state profile is found if q is sufficiently small 
(Eq.EJ). 

Startup Transients — Transients were studied by evolv- 
ing from rest using a time step of lO -5 ^ (Fig. [3]). In all 
cases shown here strongly inhomogeneous flow develops 
after the stress overshoot, leading to a sharply banded 
transient state, with a negative velocity and shear rate 
in the less viscous band. In the monotonic case the veloc- 
ity profile eventually smooths out. For a narrower stress 
plateau {e.g., (3 = 0.3, not shown) the overshoot has a 
less pronounced kink and typically a positive shear rate. 
We have found inhomogeneous transients with negative 
velocities with stress differences corresponding to a cone 
angle 9 = 0.003° in startup runs, but for 6 = 0.001° 
the amplitude of the inhomogeneous flow is reduced, and 
the velocity no longer negative, whilst for 9 = 0° the 
flow remains homogeneous. With perturbed initial con- 
ditions then the inhomogeneous transient behaviour re- 
turns. The transient for a monotonic model (/? = 0.728) 
in which spatial gradients are artificially prohibited ex- 
hibits a slower decrease after the stress overshoot than 
in the spatially resolved model (Fig. [3]); hence, inhomo- 
geneities are important when using transient data to help 
differentiate candidate constitutive models 0, . 

Large Amplitude Oscillatory Shear (LAOS) — A sinu- 
soidal spatially-averaged shear rate was applied with fre- 
quency f2 and maximum shear rate j m , and evolved from 
rest (zero stress) until any initial transients had decayed. 
We characterize the dynamics by the Weissenberg num- 
ber Wi = jrnTd and the Deborah number De = SlTd- For 
low De (frequency) we expect to recover some features of 
the steady state behavior, such as transient banding for 
Wi roughly within the non-monotonic part of the flow 
curve; while higher frequencies (high De) should pro- 
duced sharper profiles similar to the transient behavior 
in Fig. [3j since the system cannot relax before flow rever- 
sal. At the highest frequencies we expect the reversing 
dynamics to be too fast to allow an inhomogenous state. 

Fig. [4pV shows this behavior on a "Pipkin diagram" of 
Wi vs. De, for a monotonic flow curve (/3 = 0.728) in a 
slightly curved geometry. The inhomogeneous profiles in 
the banding regime (Fig. Hp) can be represented para- 
metrically in terms of shear rate and torque, (7(2/), r(y)) 
(Fig. HJ3). At the high stress regions of the cycle a por- 
tion of the sample enters the high shear rate band as 
reported experimentally [H|]. In these calculations the 
position y„ of the interface at a given strain 70/fi = 3 
varied with De as ~ (De) a , where a ~ 0.4 — 0.6, un- 
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FIG. 3: Velocity as a function of position y = q^ 1 \n(r/R\) for 
different (@,q), at shear rate log 10 7 = 1.2 (o in Fig. [2). The 
solid (red) lines indicate the steady state profile, and dashed 
lines are transient profiles at the times shown. The transient 
torque response in the spatially uniform model is shown for 
the monotonic case (j3 = 0.728) by a dashed line 



like the fixed position reported in Ref. [l5(. We suspect 
that these experiments did not attain steady state. The 
torque overshoot is typical of polymer solutions, and re- 
sembles that found in [l|| (Fig. [4j. At low frequencies 
the system has time to find a selected stress, which re- 
mains constant for part of the cycle while the shear band 
grows into the cell. At high frequencies the fluid cannot 
relax or shear band, which leads to a sinusoidal response 
and a nearly affine spatial profile. 

Step Strain — Some experiments on the relaxation after 
a step strain found a strong inhomogeneous recoil that 
developed a negative velocity gradient [13| • We illustrate 
this with a monotonic constitutive curve (/3 = 0.728), 
Fig. [3 As in Fig. 5 of [13=], an inhomogeneous velocity 
profile develops and the velocity becomes negative as the 
system recoils from the applied shear. Experimentally, 
the total displacement after recoil is of the order of a 
tenth of the gap size which is comparable with that ob- 
served here. Thus, inhomogencitics are important when 
using step strain data to help discriminate among candi- 
date constitutive models, as when using the DE damping 
function QUI . 

Summary — We have shown that behavior reminiscent 
of shear banding, as reported recently, can be reproduced 
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FIG. 4: (A) Pipkin diagram ( Wt vs De) for j3 = 0.728 and 
q — 2 x 10~ 3 , showing regions with homogeneous (o) and 
inhomogeneous profiles (•). (B) Parametric shear rate-torque 
profiles (7(3/), F(y)) (dashed: c,d) overlaying the constitutive 
curve (solid). (C) Torque evolution for different De noted on 
(A). (D) Velocity profiles for ii); profiles (a-e) are at different 
times in the cycle for De ~ 3 (ii) (o in C). 



t<0.2r d , jr d = 40 



t > 0.2r d , jT d 




FIG. 5: Relaxation of a step strain of 7 = 8, applied using 
jTd — 40 for t — 0.2rd. (Left) Velocity before the shear rate 
stops and (right) snap shots of recoil velocities. 



using the Rolic-Poly model supplemented by a term to 
accommodate spatial gradients. The RP model contains 
an unknown parameter (3, which controls the efficacy of 
convected constraint release. Even for (3 large enough 
to yield a stable (monotonic) constitutive curve, shear 
banding signatures can appear if the "stress plateau" is 
flat enough: (1) a geometry with a high stress gradient 
can induce a flow profile that could be mistaken for band- 
ing; (2) sharp banding-like profiles can appear in start-up 
transients even though the steady state is non-banded; 

(3) LAOS can trap these sharp transient profiles; and 

(4) relaxation after a large step strain can be very inho- 



mogeneous, sometimes with a negative shear rate recoil. 
Several recent experiments, particularly on polydisperse 
polymer solutions, may fall into this category [l6(. A 
wide plateau is believed to accompany very highly en- 
tangled systems Q , and the larger number of relaxation 
times are likely to render polydisperse systems intrinsi- 
cally more stable than monodisperse systems as was 
noted in recent experiments |l6[. Our results are not 
specific to the RP model; Zhou et al. recently studied a 
different two-fluid model of shear banding (with a non- 
monotonic constitutive relation) , and found qualitative 
results similar to some of ours 23 1. 
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